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Completeness of the eigenfunctions of a self-adjoint Hamiltonian, which is the basic ingredient 
of quantum mechanics, plays an important role in nuclear reaction and nuclear structure theory. 
However, until now, there was no a formal proof of the completeness of the eigenfunctions of the 
two-body Hamiltonian with the Coulomb interaction. Here we present the first formal proof of the 
completeness of the two-body Coulomb scattering wave functions for repulsive unscreened Coulomb 
potential. To prove the completeness we use the Newton's method [R. Newton, J. Math Phys., 
1, 319 (I960)]. The proof allows us to claim that the eigenfunctions of the two-body Hamiltonian 
with the potential given by the sum of the repulsive Coulomb plus short-range (nuclear) potentials 
also form a complete set. It also allows one to extend the Berggren's approach of modification of 
the complete set of the eigenfunctions by including the resonances for charged particles. We also 
demonstrate that the resonant Gamow functions with the Coulomb tail can be regularized using 
Zel'dovich's regularization method. 

PACS numbers: 03.65.Nk, 03.65.Ca, 24.30.-v, 21.60.CS 



Completeness of the eigenfunctions of a self-adjoint Hamiltonian as a part of the hnear superposition principle 
provides a powerful tool in different areas of quantum physics. Consider the set of the radial eigenfunctions S = 
Wniir), ipiik^r)}, 0<n<N,0<k<oooi the of the self-adjoined Hamiltonian H — T + Vi. Here, T is the 
kinetic energy operator, Vi{r) = (r) + Vi{r) is the sum of the short-range and centrifugal potentilas, (fniir) are the 
normalized to unity eigenfunctions of the discrete branch of the energy spectrum (bound states), n is the principal 
quantum number and I is the orbital angular momentum; ipi{k,r) are the eigenfunctions of the continious spectrum 
(scattering states) normalized to delta functions. For simplicity interacting particles are assumed to be spinless. 
Definition of the completeness: set S of the eigenfunctions of the self-adjoined Hamiltonian forms a complete set in 
the Hilbert space if any function h{r) belonging to this Hilbert space can be expanded in terms of the eigenfenctions: 
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An elegant proof of the completeness of the eigenfunctions of the two-body Hamiltonian with the interaction potential 
(r) satisfying the conditions 



dr r\V^ {r)\ < oo (6) 





and 



dr r^^+^\V^ {r)\ <oo (7) 







has been presented long ago by R. Newton P, . The first condition guarantees that a regular solution is entire 
function of k and the second one secures behavior of the eigenfunctions at /c = 0. Newton's method is based on 
consideration of the integral taken along the closed contour: 

CSO 

/(r) = j dkk j dr' h{r')G^i^\r,r']k). (8) 
(C) 

Here, G|^^(r, r';fc) is the two-body partial wave Green's function in the coordinate space. The closed integration 
contour C in Eq. ((SJ is shown in Fig. ^ It goes along the real axis in the momentum plane k bypassing the origin 
fc = along the infinitesimal contour 7, and along the semicircle R with the radius \k\B, — > 00. 

Although nobody doubts that the completeness holds also for the eigenfunctions of the Hamiltonian containing the 
unscreened Coulomb potential, a formal proof of it has not yet been presented in the literature. The reason for that 
is a quite complicated analytical behavior of the Coulomb Jost solutions fl^\k,r) For the Coulomb interaction 

fi^\k,r) {Ji \k,r)) has singularity at fc = with the cut along the negative (positive) imaginary semiaxis in the k 
plane 3]. It makes behavior of the Jost functions and Jost solutions in a vicinity of A: = quite complicated. Note 
that introducing the screened Coulomb potential or a box changes the analytical properties of the Jost solutions: the 
singularity at k = dissapears (if the potential does not support a bound state at E = 0) and the Coulomb cut along 

the imaginary semiaxis begins from k ~ —i ^ (k ~ i j^) for fl^\r) {fli ■'(r)), where R is the screening radius, as it 
takes place for short range potentials. The Coulomb analytical features can be recovered only at i? = 00. 

In this paper we present, for the first time, a formal proof of the completeness of the eigenfunctions of the radial 
two-body Hamiltonian with pure repulsive Coulomb potential at arbitrary angular orbital momentum. Our proof is 
based on the consideration of the Newton integral l] . The idea of the Newtons proof, should be applied to the pure 
Coulomb case, is as folloows: according to the Cauchy's theorem the integral over a closed contour C, Eq. JHJ, which 
does not contain any singularities of the integrand on or inside the contour is zero. Then we need to show that the 
integral over large R in he limit \k\ii — > 00 gives 7r/i(r). This result is similar to the short range case. The most 
delicate problem is to prove that the integral over small semicircle 7 disappears in the limit |fc|^ 0, where \k\j is 
the radius of the contour 7 in the momentum plane. The last step needed to prove the completeness is to show that 

the integral over the real axis can be reduced to ~n J dr'h(r') J dkipl'^^* {k, r') ^'i'^\k, r). 



II. GENERAL EQUATIONS 

From Fig. ^it is clear that the closed contour C in the complex k plane consists of the different parts: C = R + O+'y. 
Here R is the integration contour over the large semicircle, is the integration contour along the real axis with 
eliminated small interval around fc = and 7 is the integration contour over the small semicircle around k = 0. 
Correspondingly, we can split the integral in Eq. (|SJ into three parts: 

I{r)^Io{r)+I^{r)+In{r). (9) 

Here, the integral along the real axis in the fc-plane 



/o(r)=/o-(r)+/+(r). 



(10) 



Im k 
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The integrals along the contours R and 7 are given by 

00 

R 

and 

oc 

/^(r) = [dkk [ dr' h{r')Gf^'^^\r,r';k), 



correspondingly. The partial Coulomb Green's function in the configuration space can be written as 

vf\k,r^)fP^+\k,ry) 



Gf^'''^\ry-k) = 



L 



(C)(+) 

I 



(fc) 



(13) 



(14) 



(15) 



Here, fP'-^\k,r) are the singular (at r = 0) (Jost) solutions for the l-th partial wave of the radial Schrodinger 
equation with the pure Coulomb interaction 



and 



' w V ; r{l + l±iri) 

are the Coulomb Jost functions; Lpf^\k, r) is the regular (at r — Q) solution given by 

ipf\k,r) = {2ik)-^ \Lf^^^-\k)jP^+\k,r)-Lf^^+\k)fP^-\k,r) 

^ ' [T{l + l-iT]) ^ ' ' T{l + l + ir]) *" ^ ' ; 

= r'+i e''^'- ii^i(/ = 1 + ir;, 2? + 2; -2iA:r). 



(16) 
(17) 

(18) 

(19) 
(20) 



Also r/ = Zi Z2 /i/fc is the Coulomb parameter of the interacting particles 1 and 2, Zt e is the charge of particle i, /i is 
the reduced mass of particles 1 and 2. We use the systems of units such that Ti — c = 1. A physical scattering wave 
function il^l {k,r) which is normalized to delta- function is related to the regular solution (p^ {k,r) as 



L 



(C)(+) 



(fc) 



(21) 



III. INTEGRAL /(r) 



Following Newton's proof we need to show that the integral I{r) over the closed contour C, Fig. is zero. To prove 

(c 
I 



this it is enough to show that the Green's function G'['^'"'^'' (r, r' ; fc) is a regular function on the integration contour 



and inside it and that integral /(r) converges. 

According to Q the regular Coulomb solution of the Schrodinger equation (/^['^•'(fc, r) is an entire function of k in 
the complex k plane and is given by a series which absolutely and uniformly converges. The singular solution is given 

by 

/f|)(fc,r) = e^''M^_,^,,+i(-2zfcr) (22) 
where the Whittakcr function is determined in a standard way J5|]- 



I arg(-2iA:r)| < TT, (24) 
lm/c>0, -i = e-^e''^/^ ^ e-"/2 (25) 
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To select the branch of the Whittaker function in the upper half of the k plane (/to fc > 0) , where the closed contour 
C lies, we impose condition 1)24(1 on the argument of the Whittaker function. The Whittaker function is analytical 
function in k plane with the a branch point singularity at k=0 and the cut going along the negative imaginary semiaxis 
from fc = to —zoo, the blue line in Fig. Taking into account that —i = e~"e*'^/^ = e~*'^/^, we get that on the 
bank of the cut in the fourth quadrant arg k = — 7r/2, and on the left bank of the cut arg fc = 3/2 tt. At such definition 
\arg{~2ik)\ > n. Thus the Whittaker function, and, hence, the singular solution are regular functions along and 
inside the closed integration contour. We note once again that the integration contour by passes the singularity of 
the singular solution at fc = 0. Evidently that the Jost function L|^'^^j(fc) is also regular function along and inside the 
closed integration contour. The convergence of the integral follows from the behavior of the integrand in the upper 
half of the k plane. It will be shown below when we consider each integral Ir^t), Io{r) and I-y{r). Thus the Newton's 
integral 

oo 

/(r) = (j) dkk J dr' h{r')G[^\r,r';k) ^0. (26) 

(C) 

for pure repulsive Coulomb case 



IV. INTEGRAL /^(r). 

First consider the integral /_R(r) over the large semicircle R; k £ R : \k\ oo. Following Newton we split this 
integral into two parts: 



lBir)^lR,{r)+InJr), (27) 



where 



R 



Lf^^^\k) 



IrA^)^- I dkk I dr'h{r') ^' ^ l^^^'^l^ \ fcei?, |fc|->oo (28) 



and 



lR^ir) = - J dkk J dr'h{r')r^^^j^^ k € R, \k\ ^ cx^. (29) 

First we consider the integral Ir^ (r). To evaluate it, we replace the regular solution (p^ (fc, r) by its leading asymptotic 
term 

^f\k,r) "fe" (2fc)-'-iffi±^[e-*t('+i)e^('=^-^i"(2M) 

r(i + 1) 

_|_gif (/+l)g-i(fer-j)ln(2/£r))j^ ^30) 

Similarly the leading asymptotic term of the singular solution is given by 

//^)(+)(fc,r) (31) 
Finally the asymptotic behavior of the Jost function is given by 

4S,ik) (2fc)-e-f (32) 

Taking into account the expression for the partial- wave Coulomb Green's function, Eq. ((15(1 . we get its leading 
asymptotic term 

G'P(+)(r, r-'; k) = i{2k)-^ [(-1)' e''^'(''+''') - e^'^^'^''')] . (33) 
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Substituting it into the integral (|29|l gives 



Ij^^{r)=i dkk— / dr'h{r') (^-^ly e'Hr+r') _ ^rHr-r') 



(34) 



We evaluate each term separately using integration by parts. The first term containing e*'^(''+'' ^ tends to zero for 
oo. The second term containing e*'^^'' reduces to 



Similarly we get that 



(35) 



Hence the total integral /^(r) reduces to 



(36) 



(37) 



Thus we arrived exactly at the same result as for the short range interaction 2]. It is not surprising because the 
leading asymptotic term of the partial Coulomb Green's function has no trace of the Coulomb interaction and behaves 
like a leading asymptotic term of the free Green's function. 



V. INTEGRAL /o(r) 



The integral over the real axis Io{r) is given by the sum of two integrals, /g \r) and l'^\r), Eqs Hll|l and 
(|12|l . We split into two parts as we did for lR,{r): 



(38) 



(39) 



Here 



±oo I \ ) 



(40) 



and 



itoo oo 



/±^(r) = -lim dkk dr'h{r') 



Lf^^+\k) 



(41) 



Taking into account that ip\ (A;, r') — ipj^ (~^, ^0 we can rewrite the integral /o< (^) 



IoAr)^-J dr'h{r') 





YunlAkkkJ^Hk fr'^\-k,r) 



(42) 



However 



fr'^\-ky) fr~\ky) 



i-k) 



\k) 



(43) 
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Then the integral (|42|l reduces to 



/ojr) =- / dr'h{r') 



= - J dr' h{r- 



Taking into account Eq. (|18|l we get for Iq^ (r) 

r 

/oJr) = -2i Jdr'hir') 



Similarly we get for /q^ (r): 

oo 

/o,(r) = -2i Jdr'hir') 

r 

Then the integral /o(^) becomes 

CSO 

/o(r) = -2i [ dr'h{r') 



hm / dfcfc^f)(fc,r')^^^ ^ ' ' \cJ) fcH+) ' ^ 



(44) 
(45) 



lim / dfc k 



,^Pik,r')vP{k,r) 



(46) 



,ipf\k,r')v\^>ik,r) 



lim / dfc fc 

e-+0 



Lf(fc) 



(47) 



lim / dfc fc^ 

£^0 



ipf\k,r')^f\k,r) 



(48) 



Let us consider now the limit e — > 0. As we have mentioned the regular solution (pi {k,r) is the entire function of 
fc in the finite complex fc plane 4:j. It is given by a series which is absolutely and uniformly converges in the finite fc 
plane (JfjL e. \ip'"i'^\k,r)\ < A{r) in the finite complex fc plane including fc = 0. Now let us consider l/|L^''(fc)p for 
fc^+0[l2: 



l/|Lf (fc)p- 



1 



-77 rj 



2'Kr] 



[(2? + l)!]2 e-^^-e-"" 
Thus we can take the limit e +0 in Eq. H48(l getting 



-{2kr\{{f + r, 



2 , 2n fe^+o e' 



-2 IT rj 



(49) 



/o(r) = -2i J dr'h{r') 



dfcfc 



,^f)(fc,r')^r)(fc,r) 



(50) 



Again, for the integral along the real axis we got the same results as for the short range interaction. It is not surprisng 
because eventually we used the symmetry properties of the regular solution under fc — > — fc, which is the same as for 
the short range interaction, and limiting behavior of the integrand for fc — > +0. 



VI. INTEGRAL /^(r). 



The integral over the small semicircle around fc = is the most difficult part of the proof: the point fc = is a 
regular point for well behaved short range potentials and singular point for the Coulomb potential. We cannot avoid 
this integral by just writting /(r) = Ir{t) + lim /□(t', e) and taking the limit e ^ because the integration contour C 

e— +0 

in this case is not a closed contour and we cannot use Cauchy theorem to evaluate I{r). We split again the integral 
/^(r) into two parts: 

/^(r) = /^Jr)+/^^(r). (51) 
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Here 



and 



ir)=- Idkk I dr' hir') ^^''^l^^y^'^ (52) 



I^^ {r) = - j dkk j dr' h{r' 
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We will prove that — > as the radius \k\^ of the semicircle 7 goes to zero. We start from consideration of 

Xy^(r). The regular solution ipi^\k,r) does not generate any problems for A; — > for Imk > because it is entire 
function of k in the finite complex k plane. The singular solution can be written in the form: 

00 

/f )W(fc,r) = e--'/^ r(; + l + i^) "''^ (-2^^^"' / d*e-*t^' (1 - ^)--+'. (54) 



Letting z = 1 — 2 ikr/t we get 

oc 

fi^^^^\k,r)=e'"^^^——^ —e''''^{-2ikr)-' I dte-H^^ z'^^^+K (55) 





We can write down 



(56) 



Here, 



z=\z\e''', 6> = - arctan "^^^^ , (57) 

' ' t + 2k2r ^ ' 

A;=|fc|e'"^, ^ = arctan ^ = 7r/2 — arctan -i. (58) 

k\ ki 

It is evident that — 7r/2 < 9 < 7r/2 and < </> < tt because —e<ki = Rek < e and k2 = Imk > for fc e 7. 
Taking into account that 



,/c2r , ^kfr"^ Jk\r , Jk\'^r'^ 



(l + 2-)^+4^<yi + 4^+4 

?Ei:>i 

and 1^ sin0 > for ^2 > we get < 1. Besides 

. a , a , 2k\T , . , /ci , 

U — cos© = — — arctani ; — ) sm( arctan — ) 

k ^ k ^t + 2k2r' ^ k2 



< 1 + tl^ > 1 (59) 



= --arctan( — ' < 0, kGj. (60) 

k H + 2k2r' ^/^TM - ' ' ^ ^ 



Then 



Correspondingly for fj'^^^~^\k,r) we get 

00 

|/r^+^(fc,r)| < le^'^ ^^i^l^.^^ (-2ifcr)-'| y dte'^^ 



00 

^l^^'^' ^Z + l + ir?) (-^^fe^r'l Jdte-H^{t + 2\k\ry. (62) 
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Note that in the integrand we majorized |z|' by [{t + 2 |fc| r)/t]'following Eq. (|59|l . We also took into account that for 
k G J {k2 > 0) e*'^''! < 1. It is evident that the integral over t converges and can be uniformly majorized in the finite 
k plane. Let 



J{\k\r) = / dte-'f {t + 2\k\r) 



(63) 



Let s be a semicircle in the upper half k plane with the radius |fc|s > |fc|^. Then for all fc € s J{\k\ r) < J{\k\s r) and 
'^^^ majorized as follows: 



p(C)(+)n. '''1^° I.T,,/2 L 



Taking into account Eq. (|17|l we get 



r(2L + 2) 



(64) 



(65) 



Then it is straigtforward to see from Eq. (|52(l that 

|/^,(r)| = | Jdkk Jdr'Hr') ^^'"^^ ^'"'^icxiT^^^ -^(1^1-^) 1^1 I J dr' h{r') A{r')\ ^ 0, (66) 



7 

-(C) 



where A{r) > |e[ {k,r)\ is the majorizing function for fees. We also took into account that 7 is | / dfc fc | = 2 



where 7 is the semircle. Similarly we prove that 



^f)(fc,r)/f)W(fc,r'), 



IA>(^)I = \ dkk / dr'h{r') 



Lf^^+\k) 



< 2 



1 



r(2L + 2) 



A{r) \k\ I / dr' h{r') J(|fc|, r')\ 0. (67) 



Note that assumption that h{r) £ is enough to provide the convergence of the radial integral in Eq. (|67|l . 
Thus we proved that |/^(r)| < ||/^<(r)| + |/^>(r)| '''J^" 0. 



VII. COMPLETENESS OF THE COULOMB SCATTERING WAVE FUNCTIONS. 

Now we can return to Eq. Replacing /(?') and I.y{r) by zero, /^(r) by Eq. (p?7jl and loir) by Eq. (pjO]) we get 



h{r) = - y dr'h{r') 




dkk' ^\ky)^f\k,r) 



LfHk) 



(68) 



From this equation we may conclude that the Coulomb scattering wave functions for repulsive Coulomb interaction 
satisfy the completeness relationship: 



00 

S{r' -r) = - J dkk^ 



^P{k,r')^f\k,r) 



Lf\k) 



(69) 



Introducing the spectral function p{E) we can rewrite the completeness relationship (|69f) in the form: 

00 

5{r' -r)= J dp{E) ipf^ {k, r') ifP {k, r), (70) 
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where 



dp j 2_A£fe 
= < 7r 

dE I n 



E>0, 
E <0. 



(71) 



Using the physical scattering wave function H21() we can rewrite completeness relationship ()69|l in the standard form 



S{r'-r) = - I dk^P*{k,r')iP\'^'>{k,r). 



(72) 



Thus for the repulsive Coulomb interaction we got the same result as Newton Q, 0] for the short range interactions: 
the system of the physical scattering wave functions S ~ {"0; {k,r)} forms a complete set. 



VIII. BERGGREN'S METHOD AND NORMALIZATION OF THE RESONANCE GAMOW 

FUNCTIONS FOR CHARGED PARTICLES 

The completeness of the eigenfunctions of the two-body Hamiltonian with the short range interaction and the 
formal proof of the completeness of the Coulomb scattering wave functions presented here allows us to claim that the 
eigenfunctions of the Hamiltonian for the repulsive Coulomb + short range (nuclear) potential form a complete set 

^ i \Lf\k) 

Here e„; is the normalized to unity bound state wave function with the principal quantum numcr n in the partial wave 
I. In Berggren's method 6] the complete set of the bound states and continuum states for real positive energies E > Q 
can be redefined by including the resonant states. A new complete system consists of the discrete states, bound and 
resonant, and continuum states, scattering states for real positive and complex energies. To normalize the resonant 
states Berggren intdoduced the dual basis and used Zel'dovich's regularization procedure 0. However, Zel'dovich's 
regularization procedure has been shown to work only for the resonances for the short range potentials. We will show 
here that Zel'dovich's regularization works also for the resonant Gamow states with the Coulomb tail. 

Let ip^iif) stands for the Gamow wave function describing the resonant state in a system of two charged particles 
interacting via the sum of the repulsive Coulomb + nuclear potential. The dual Gamow function (f^^ (r) satisfies the 
condition: ip^j' (r) = (p^^ (r) . Asymptotic behavior of the Gamow function is given by 

(74) 

where bni is the amplitude of the tail of the resonant Gamow function, k^ — ki + ik2 is the momentum of the 
resonance, fci = Rek > and ^2 = Imk < 0; -q^ = Zi Z2 n/k^ is the Coulomb parameter for the resonant state. 
The radial factor ^ in the denominator appears due to the Coulomb barrier. Zeldovich's normalization condition 
of the Gamow function is given by 

00 

iV = lim 1 dr ^« (r) (r) e-^^\ (75) 


The introduction of the regularization factor exp(— /3r^) is neccessary to provide the convergence of the normalization 
integral for r ^ 00 which otherwise diverges because exp(z k^ r) = exp(« k\r — k^r). We split the integral ((75|) into 
two parts: 

A 00 p, 

N = hm [ / dr ^« (r) (r) + 6„z / dr e"'^ 





A 



dr (r) cp^t (r) + lim / dr e"^ . (76) 



A 



11 



A is assumed to be large enough to approximate at r > A the Gamow function by its leading asymptotic term H74|l . 
Evidently, the integral over a finite interval r < A converges and we can take limit P —^ in this integral. We need 
to prove that the second integral can also be determined in the limit /3 — > 0. Taking into account that 



^""^ 1^ = (fcO^TT^^^^ -^k,) = X-^S, (77) 

where a = Zi Z2 ^, and A > and (5 < 0. Hence 

00 

7Ve.t = bni hm / dr 2 '^^ '^^ ' * e-'' '■^ (78) 

A 

For neutral particles there is only one exponentially diverging factor exp(— 2 fc2 (for r — > oo) in the integrand. 
Zel'dovich's regularization method was applied to regularize such integrals. This procedure works for ki > The 
presence of the oscillating exponential exp(j 2 ki r) is crucial for Zel'dovich's regularization. However for charged 
particles due to the Coulomb barrier at complex momentum we have an additional diverging factor in the integrand, 
7--2<5^ We will show now that Zel'dovich method works for the charged particles also. We rewrite Eq. (|78() as the sum 
of two terms: 

A CO 

N,,t = N^ + N2^ -bra lim / dr ^ fc- . ^-2 . ^-/j ^ ^^^^ / dr ^ ^ r"^ ^ n'^ . (79) 



The first term in this equation converges and we can take limit /3 ^ in this term. It is enough to consider the 
second term only. The second integral can be taken analytically using Eqs. (3.462(1)) and (9.246) [|: 

OO 

lim /dre'2fcr^-2.^ -/3r^ ^ lim (2/3)^r(-2i77+l)e-^^-D2»„-i(-^) 



= r(-2i?7+l)(-i2fc)2'"'-i =J{k). (80) 

Here, D^{x) is the parabolic cylinder function. For Imk < and Re{—2iri) > the integral on the left hand-side of 
this equation determined only for /3 > 0. Meantime function J{k) is an analytical function of k in the whole complex 
k plane except for the singular points at 2 i 77 = n, where n — 1,2..., and fc = 0. Hence function J{k) can be considered 

oc 

as an analytical continuation of the integral J ^7. e« 2 fcr ^,-2 ir; -^^^^ region Imk < 0. Thus we have shown that the 



Zel'dovich's regularization can be used to determine the norm of the Gamow functions for charged particles. 



IX. SUMMARY 



We presented a formal proof of the completeness of the eigenfunctions of the two-body Hamiltonian with repulsive 
Coulomb potential using the Newton's integral P, 01 • The most delicate point was to investigate the behavior of 
the integrand near the singular point A; = and to prove that the contribution from a small semicircle 7 around 
k = goes to zero as its radius |fc|^ —f 0. The presented proof allows us to claim that a system of the eigenfunctions 
S — {(pni{r), tpi{k,r)}, 0<n<N,0<k< 00, of the of the self-adjoined Hamiltonian with the potential given by 
the sum of the Coulomb and nuclear interactions is also complete. 

Recently the Berggren's method 6] has been used in the so-called Gamow shell model 8]. The Berggren's method 
can also be used for the Green's function spectral decomposition in nuclear reaction theory with charged particles 
. Our proof validates application of the Berggren's method for charged particles. We have also demonstrated 
that the Gamow resonant states with the Coulomb tail can be normalized using Zel'dovich's regularization method 
0. It is a crucial point for application of the Berggren's technique. Note that Zel'dovich's regularization method is 
not unique and other regularization techniques can be used. 
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